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INTRODUCTION 


Ocean  Thermal  Energy  Conversion  (OTEC)  is  the  concept  for  a  system 
which  will  extract  energy  from  the  ocean  by  taking  advantage  of  the 
sizable  temperature  difference  between  surface  and  bottom  water  in  a 
tropical  ocean^y  .  A  fundamental  part  of  an  OTEC  system  is  a  large  sec¬ 
tion  of  heat  exchangers  through  which  the  energy  transfer  is  made.  It 
is  crucial  to  keep  the  surface  of  these  heat  exchangers  clean  because 
fouling  by  organisms  in  sea  water  can  cause  deterioration  of  the  heat 
transfer  effectiveness.  The  Naval  Coastal  Systems  Center  (NCSC)  inves¬ 
tigated  cleaning  techniques  that  would  effectively  prevent  fouling  in 
metal  tubes  over  a  long  period  of  time. 


The  physical  structure  and  cleaning  systems  of  the  OTEC, assembly  at 
Panama  City,  Florida  have  been  described  in  several  papers'^  ^ V  and 
will  be  dealt  with  only  briefly  here.  This  report  will  describe  the  com¬ 
puter  driven  control  and  data  analysis  system  used  in  conjunction  with 
the  mechanical  system.  This  will  include  explanations  of  the  physical 
computer  setup  and  its  interaction  with  hardware  components  of  the  sys¬ 
tem,  the  software  methods  of  sampling  and  data  management,  and  data 
analysis  techniques.  A  brief  discussion  of  the  underlying  theory  of 
heat  transfer  is  also  included. 


(  ^Griffin,  0.  M. ,  OTEC,  Power  from  Thernut\C,radients ,  Sea  Technology, 
pp.  11-15,  August  1977. 

^2'Frit8ch,  A.,  Adamson,  W. ,  and  Casteili,  V.,  An  Evaluation  of  Mechani¬ 
cal  Cleaning  Methods  for  Removal  of  Soft  Fouling  from  Heat  Exchange r 
Tubes  in  OTEC  Power  Plants,  Proceedings  of  the  Ocean  Thermal  Energy 
Conversion  (OTEC)  Biofouling  and  Corrosion  Symposium,  Seattle, 
Washington,  1978. 


(^Braswell,  J.  A.,  Lott,  D.  F.,  and  Hedlicka,  S.  M.,  Preliminary 

Evaluation  of  Flow=  Driven  Brushes  for  Removal  of  Soft  Bio fouling  from 
Heat  Exchanger  Tubes  in  OTEC  Power  Plants,  Ocean  Thermal  Energy  Con¬ 
version  (OTEC)  Workshop,  Washington,  D.  C.,  January  1979. 


^Lott,  D«  F.,  and  Tuovila,  S.  M.,  Biofouling  Countermeasures  -  Status 
of  Two  Mechanical  Systems  Systems  and  Chlorination,  Sixth  Ocean  Thermal 
Energy  Conversion  (OTEC)  Symposium,  Washington,  D.  C.,  June  1979. 
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TEST  SITE  CONFIGURATION 


The  test  site  at  NCSC  is  located  on  the  shores  of  St.  Andrew  Bay, 
an  estuary  of  the  Gulf  of  Mexico,  at  Panama  City,  Florida.  The  OTEC 
piping  assembly  and  computer  system  was  located  on  a  pier,  with  the 
seawater  being  pumped  to  test  sections  from  a  depth  of  approximately 
7  feet^3).  Three  self-priming  centrifugal  pumps  were  available  with 
one  or  two  being  used  to  pump  seawater  from  the  bay  into  a  header. 
Twelve  tubes  were  fed  from  the  header;  six  were  of  titanium  and  six  of 
5052  aluminum  so  that  the  cleaning  systems  could  be  tested  on  two  dif¬ 
ferent  metals.  Water  flowed  from  each  tube  back  into  the  bay,  with  the 
flowrate  manually  adjusted  with  a  valve  at  the  outlet. 


CLEANING  SYSTEMS 


The  tubing  system  contained  four  control  tubes,  two  tubes  in  a  flow 
driven  brush  cleaning  system,  two  tubes  in  a  reciprocating  sponge  ball 
cleaning  system,  and  four  tubes  in  a  chlorination  system.  All  cleaning 
systems  worked  automatically  and  were  controlled  by  timers. 

One  pair  of  control  tubes  (one  aluminum  and  one  titanium)  were 
allowed  to  foul  freely  to  provide  data  on  the  amount  and  kind  of  fouling 
present.  The  other  pair  of  control  tubes  were  cleaned  daily  and  served 
as  an  internal  check  of  the  data  gathering  system. 

The  flow-driven  brush  cleaning  system  consisted  of  one  aluminum 
and  one  titanium  tube  through  which  a  nylon  brush  was  passed.  At  each 
end  of  the  tube  was  a  nylon  basket  in  which  the  brush  was  trapped. 
Periodically,  water  flow  in  the  tube  was  reversed  to  force  the  brush  to 
the  other  end  of  the  tube.  The  brush  diameter  was  slightly  larger  than 
the  tube's  inner  diameter  and  would  shear  the  fouling  material  off  the 
tube  walls. 

The  recirculating  sponge  ball  cleaning  system  consisted  of  one 
aluminum  and  one  titanium  tube  in  which  a  sponge  ball  of  a  slightly  lar¬ 
ger  diameter  than  the  tube  was  circulated.  The  ball  was  held  in  a  valve 
and  released  periodically  when  the  valve  was  turned.  When  the  ball  had 
gone  past  an  optical  sensor,  the  valve  was  turned  back  to  its  original 
position  and  the  ball  was  diverted  into  a  catcher  which  placed  it  back 
into  the  circulating  system. 


ibid. 
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The  control,  brush,  and  sponge  ball  systems  have  been  described 
by  Lott  and  Tuovila^4^.  Tests  with  the  chlorination  system  have  nor 
yet  begun.  This  system  consists  of  two  aluminum  tubes  and  two  titanium 
tubes.  A  chlorine  generator  releases  a  known  concentration  of  chlorine 
into  the  system;  chlorine  release  may  be  pulsed  at  regular  intervals  or 
may  be  continuous  so  that  the  concentration  of  chlorine  within  the  tubes 
remains  constant. 


COMPUTER  CONTROL  SYSTEM 


A  DEC  PDP-11/34  computer  was  used  to  provide  system  control  and 
data  analysis.  Ideally,  this  system  should  run  continuously  without 
much  operator  assistance.  The  main  areas  of  computer  control  are  pump 
selection,  tube  selection,  and  heater  control;  if  necessary,  the  entire 
test  could  be  halted  by  the  computer. 

Pump  performance  was  monitored  throughout  a  test.  If  a  pump 
failed  to  provide  uniform  flow,  an  error  message  was  printed  to  alert 
the  operator  to  a  problem.  If  the  bad  pump  performance  continued,  the 
computer  was  programmed  to  automatically  shut  down  the  malfunctioning 
pump  and  simultaneously  activate  another  pump  through  a  relay  panel. 

Heaters  in  the  OTEC  system  consisted  of  copper  heater  cylinders 
which  were  press  fitted  to  the  outside  of  the  heat  exchanger  tubes. 
Voltage  to  the  heaters  was  controlled  by  a  Variac  panel,  which  must  be 
set  manually  and  read  into  the  computer  via  a  DEC  ADll-K  analog  to  digi¬ 
tal  converter.  The  computer  was  programmed  to  be  responsible  for  turn¬ 
ing  individual  heaters  on  and  off,  monitoring  rates  of  heat  up  and  cool 
down  of  sea  water  in  the  tubes,  and  setting  heater  voltage  limits.  The 
timing  schedule  for  heat  up  and  cool  down,  and  the  allowable  voltage 
limits  were  pre-set  by  the  computer  but  these  default  conditions  could 
be  altered  by  the  operator  at  the  beginning  of  a  test.  Once  a  test  had 
begun,  the  computer  would  continuously  turn  the  heaters  on  and  off  ac¬ 
cording  to  the  prescribed  cycle,  monitor  heater  performance,  and  peri¬ 
odically  sample  and  print  voltages.  If  a  heater  built  up  more  voltage 
than  a  prescribed  upper  limit,  that  heater  was  turned  off  by  the  com¬ 
puter  and  the  corresponding  tube  was  "dropped"  from  the  test.  This  was 
to  prevent  damage  to  heat  transfer  units. 


'4Hbid. 
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Sonic  flow  meters  were  used  to  measure  flow  velocity  and  read  into 
the  computer  via  the  A  to  D  converter.  Flow  had  to  be  set  manually  via 
a  PVC  valve  at  the  output  of  each  tube,  but  the  computer  monitors  flow 
rate  throughout  a  test.  Maintaining  a  constant  water  flow  was  crucial 
to  correctly  determining  the  heat  transfer  characteristics  of  a  tube, 
so  any  sudden  change  in  flow  would  trigger  an  error  message  to  alert 
the  operator.  Intermittent  or  sluggish  water  flow  would  cause  the 
heaters  to  heat  the  tube  excessively  or  erratically,  and  that  tube  would 
be  dropped  from  the  test. 


DATA  SAMPLING 


Software  sampling  programs,  in  FORTRAN  IV  and  MACRO  assembly  lan¬ 
guage,  were  initially  planned  and  written  by  personnel  at  the  David  W. 
Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC)  at  Annapolis, 
Maryland.  These  routines  allowed  the  computer  to  read  signal  output  by 
hardware  components  of  the  OTEC  systems.  An  analysis  cycle  consisted 
of  two  parts:  (1)  heat  up,  when  the  heaters  are  on,  and  (2)  cool  down 
after  the  heaters  have  been  shut  off.  Samples  were  taken  during  the 
cool  down  cycle  to  determine  the  exponential  heat  decay  curve  for  each 
tube. 


The  general  data  analysis  scheme  is  outlined  in  Figure  1.  Four 
quantities  were  sampled  for  each  tube:  (1)  time,  (2)  heater  voltage, 

(3)  water  temperature  in  degrees  Celsius,  and  (4)  flow  rate  in  feet  per 
second.  Aluminum  tubes  were  sampled  successively  with  heater  voltage 
and  water  temperature  sampled  every  3  seconds,  and  flow  rate  measured 
every  9  seconds.  For  titanium  tubes,  the  sampling  rates  were  twice 
those  for  aluminum  tubes  because  the  cooling  curve  for  a  titanium  tube 
has  a  longer  time  constant  than  that  of  an  aluminum  tube.  For  clean 
tubes,  a  time  constant  of  50-55  seconds  is  typical  of  the  cooling  curve 
for  a  titanium  tube  and  30-35  seconds  for  aluminum.  Sampled  values  were 
displayed  on  a  CRT,  which  was  updated  every  3  seconds,  and  printed  on  a 
teletype  every  10  minutes.  Time  was  provided  continuously  by  a  battery 
driven  quartz  clock  inside  the  computer  and  provided  for  automatic  re¬ 
starting  of  an  OTEC  test  following  any  interruption  of  power  to  the  com¬ 
puter.  Sampled  values  were  stored  in  an  array  for  the  duration  of 
sampling  cool  down  time  and  were  printed  on  a  teletype  during  the  next 
heat  up  cycle;  and  were  used  to  provide  an  analysis  of  that  tube.  Keat 
up  time  was  set  at  18  minutes  with  a  cool  down  time  of  12  minutes  yielding 
48  analyses  per  tube  each  day. 
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FIGURE  1.  OTEC  DATA  SAMPLING  AND  ANALYSIS 
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DATA  ANALYSIS 


Data  analysis  techniques  were  based  on  those  developed  at  Carnegie- 
Mellon  University^1’)  and  were  intended  to  determine  changes  in  the  foul¬ 
ing  resistance  based  on  the  transfer  of  heat  from  the  tube  walls  to  the 
sea  water  flowing  inside. 

In  determining  the  heat  transfer  coefficient  in  a  tube,  it  was 
assumed  that  all  resistance  in  the  system  was  due  to  buildup  of  a  foul¬ 
ing  layer  at  the  interface  of  the  pipe  interior  and  the  flowing  sea 
water.  The  resistance  at  other  sites  within  the  unit,  such  as  that 
present  in  the  copper  block  and  at  the  interface  between  pipe  and  copper 
block,  were  assumed  to  contribute  only  negligibly  to  the  total  resistance 
in  the  system.  ^ 


As  a  test  progressed,  the  amount  of  fouling  present  was  estimated 
from  the  changes  in  resistance  from  a  baseline  value  determined  from  a 
Wilson  Plot<6>  done  on  data  taken  from  each  tube  when  it  was  clean. 

This  plot  is  essentially  a  linear  regression  relating  the  inverse  of 
flow  velocity  (l/v0*8)  with  fouling  resistance  (1/h) ,  where  h  is  a  mea¬ 
sure  of  heat  transfer.  A  sample  Wilson  Plot  is  shown  in  Figure  2. 

Slope  was  determined  experimentally  to  be  approximately  3.4&  x  10“3  and 
the  intercept  at  l/V^.8  =  0  should  ideally  be  zero.  The  intercept  nor¬ 
mally  takes  a  small  positive  value  (approximately  1.0  x  10~^)  which 
represents  the  total  amount  of  resistance,  in  the  system  which  is  not 
due  to  biological  fouling  (assuming  the  tube  has  been  properly  cleaned) . 
As  fouling  increases,  the  Wilson  line  slope  should  remain  unchanged 
while  the  intercept  increases,  providing  an  estimate  of  how  much  the 
fouling  resistance  has  increased  over  the  baseline  value. 


During  a  test,  sea  water  flowing  through  tubes  in  the  heat  trans¬ 
fer  units  was  heated  until  the  temperature  of  the  tube  walls  stabilized 
at  a  temperature  slightly  above  that  of  the  water  flowing  inside.  Then 
the  heaters  were  turned  off  and  the  tubes  cooled  by  the  colder  sea  water 
flowing  through  them,  assuming  a  constant  flow  velocity,  the  voltage 
decay  of  the  cooling  curve  can  be  defined  by  the  relation: 

V(t)  =  V  e~t/T  (1) 

o 


^^Fetkovitch,  J.  G. ,  Fette,  C.  W. ,  Findley,  R.  W. ,  Grannemann,  G.  R. , 
Mahalingham,  L.  M.,  Meier,  D.  L.,  and  Runco,  P.  D.,  A  System  for 
Measuring  the  Effect  of  Fouling  on  Heat  Transfer  Under  Simulated  OTEC 
Conditions ,  Report  COO-4041-10,  Camegie-Mellon  University,  December 
1976. 

Wilson,  E.  E.,  A  Basis  for  Rational  Design  of  Heat  Transfer  Appara¬ 
tus,  American  Society  of  Mechanical  Engineers  Transactions,  Vo],  37, 
ASME  Paper  #1477,  1915. 

(Text  Continued  on  Page  8) 
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FIGURE  2.  WILSON  PLOT 
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where  t  is  time;  V  Is  the  steady  state  voltage  achieved  after  a  ttibc 
has  been  heated;  i°is  the  time  constant  of  the  cooling  curve;  and  V(t) 
is  the  voltage  at  time  t.  A  plot  of  an  exponential  decay  curve  Is  shown 
in  Figure  3.  As  fouling  increases  in  a  tube,  the  above  relation  will 
continue  to  hold,  with  increasingly  larger  values  of  t. 

In  actual  practice,  the  raw  cooling  curve  data  arc  modified  before 
a  time  constant  is  determined.  The  first  few  seconds  of  the  curve  are 
discarded  to  ensure  that  the  heater  was  turned  off  before  data  were  taken. 
The  actual  time  was  deduced  from  experience  with  the  physical  system  and 
set  at  15  seconds  for  aluminum  tubes  and  30  seconds  for  titanium  tubes 
for  our  system.  Secondly,  the  cooling  curve  was  normalized  to  zero  by 
subtracting  the  minimum  voltage  achieved  by  the  heater;  this  minimum  was 
estimated  by  averaging  the  values  of  the  last  several  voltages  at  the  end 
of  the  cooling  curve  where  the  voltage  had  stabilized  to  a  constant 
value.  Finally,  the  curve  was  weighted  to  minimize  bia9  caused  by  equip¬ 
ment  noise,  flow  irregularities,  and  other  random  system  perturbations. 


The  value  of  this  background  noise  bias  was  assumed  to  be  constant 
and  was  estimated  by  the  standard  deviation  in  the  points  vised  to  deter¬ 
mine  the  minimum  heater  voltage.  Since  low  voltage  values  in  this  cool¬ 
ing  curve  will  be  affected  by  the  bias  more  than  greater  voltage  values, 
the  weighting  function,  based  on  the  determined  system  bias,  weights 
higher  voltages  more  than  lower  ones,  using  the  relation: 


w(t) 


(V  (t))2 
n 

s 


(2) 


where  t  is  the  time,  Vn(t)  is  the  normalized  voltage  at  time  t.  s  is  the 
standard  deviation  of  the  "zero"  voltage,  and  w(t)  ts  the  weight  applied 
to  the  voltage  at  time  t. 

Rather  than  performing  an  exponential  regression  on  the  cooling  curve 
data,  the  logarithm,  base  e,  of  each  data  point  was  taken  and  a  linear 
equation  fit  to  the  resulting  line.  Weighting  of  the  curve  was  done  in 
the  linear  regression  routine  LINFIT^7^.  Taking  the  logarithm  of  both 
sides  of  Equation  (1)  yielded  the  relation: 

in  V(t)  -  in  V  -  -  . 

O  T 


v  '  Bev Ington,  F,  R.  ,  Data  Reduction  and  Error  Analysis  for  the  Physica  l 
Sciences,  pp.  104-103,  1969. 


(Text  Continued  on  Page  10) 
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FIGURE  3.  COOLING  CURVE  PLOT 
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Linear  regression  then  determined  the  values  of  the  unknowns  V  (starting 
or  maximum  voltage)  and  x  (time  constant  of  the  cooling  curve).  A  corre¬ 
lation  coefficient,  R,  and  root  mean  square  error,  RMS,  were  also  calcu¬ 
lated  for  the  linear  fit. 

Linear  regression  was  performed  twice  for  each  cooling  curve.  The 
fit  was  done  on  the  entire  cooling  curve  (except  for  the  points  at  the 
beginning  which  were  previously  deleted)  and  a  preliminary  estimate  of 
the  time  constant  (t)  determined.  Then  one  time  constant's  worth  of  data 
from  the  beginning  of  the  curve  was  fit  and  a  final  value  of  x  was  deter¬ 
mined.  A  more  accurate  fit  was  gained  by  not  including  the  "tail"  of 
the  cooling  curve  where  the  voltage  was  approaching  its  constant  minimum 
value. 

Once  the  value  of  x  has  been  determined,  the  heat  transfer  coeffi¬ 
cient  can  be  calculated  with  the  equation: 


2  3 

£n  h  =  A  +  B(£n  x)  +  C(£n  x)  +  D(£n  x) 


(A) 


where  A,  B,  C,  and  D  are  constants  based  on  physical  characteristics  of 
the  tube  (dimension,  heat  capacity,  thermal  conductivity,  etc.).  Table  1 
gives  a  description  of  tube  parameters  with  actual  values  for  each  tube 
being  listed  in  Table  2. 


This  value  of  h  had  to  be  corrected  because  its  calculation  was 
based  on  the  assumption  that  all  the  heat  lost  from  the  tube  walls 
during  cool  down  cycles  was  transferred  to  the  water  flowing  through 
the  tube.  Two  other  sources  of  heat  loss  were  considered:  (1)  into 
the  air  and  (2)  axial  loss  along  the  tube  wall.  Heat  loss  into  the 
air  was  allowed  for  with  the  equation:^ 8>ts)(io) 


^ Sparrow,  E.  M.,  Hallman,  and  Siegel,  T.  M. ,  Turbulent  Heat  Transfer 
in  the  Thermal  Entrance  Region  of  a  Pipe  with  Uniform  Heat  Flux, 
Appl.  Sc.i. ,  Res.,  Section  A.,  Vol.  7,  1955,  pp.  37-52. 

Hartnett,  J.  P.,  Experimental  Determination  of  the  Thermal  Entrance 
Length  for  the  Flow  of  Water  and  Oil  in  Circular  Pipes ,  Transactions 
of  the  ASME,  Vol.  77,  1955,  pp.  1211-1220. 

^10^Allen,  R.  W. ,  and  Eckert,  E.  R.  G.,  Friction  and  Heat  Transfer 
Measurements  to  Turbulent  Pipe  Flow  of  Water  (Pr  =  7  and  8)  at 
Uniform  Wall  Heat  Flux,  Trans.  ASME,  Journal  of  Heat  Transfer,  Vol. 
86,  Ser .  C,  No.  3,  August  1964,  pp.  301-310. 


(Text  Continued -on- Pago  1.3) 
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TABLE  1 

DESCRIPTION  OF  TUBE  PARAMETERS 


Parameter 

No. 

Description 

1 

A:  First  constant  for  calculation  of  h 

2 

B:  Second  constant  for  calculation  of  h 

3 

C:  Third  constant  for  calculation  of  h 

4 

D:  Fourth  constant  for  calculation  of  h 

5 

RADI:  Inside  isadius  of  tube  (in) 

6 

RAD2:  1/2  Machined  O.D.  of  tube  or  copper  heater  I.D.  (in) 

7 

RAD3:  1/2  O.D.  of  copper  heater  section  (in) 

8 

AKCYL:  Thermal  conductivity  of  copper  heater  (B/hr  ft  °F) 

9 

AKTUB:  Thermal  conductivity  of  tube  (B/hr  ft  °F) 

10 

HAIR:  Air  heat  transfer  coefficient  for  heat  loss  to  air 

11 

HINTR:  1/h  intercept  on  Wilson  Plot  (different  for  each 

tube) 

12 

ALENGTH:  Length  of  heater  section  on  tube  (in) 

13 

TTHWLL:  Wall  thickness  of  thin  tube  section  near 

heater  (in) 

14 

VNOM:  Nominal  water  velocity  to  which  results  are 

corrected  (fps) 

15 

HRSLPE:  Slope  from  Wilson  Plot  (different  for  each  tube) 

16 

THINT:  Thermistor  calibration  curve  intercept  (°C) 

17 

THSLP:  Thermistor  calibration  curve  slope 

18 

FMINT:  Flow  meter  calibration  curve  intercept  (gpm) 

19 

FMSLP  Flow  meter  calibration  curve  slope 

20 

IWNST:  Points  after  cooling  curve  starts  for  data 

collection  to  start 

21 

IWN:  Number  of  t's  for  second  data  fit 

1-19  are  real  numbers 
20-21  are  integers 
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TAP.LE  2 

VALUES  OF  THE  PARAMETERS 


Parameter  _ Tube  No. 


No. 

1 

2 

3 

4 

5 

0 

I 

U.  577000 

49.278000 

13.577000 

49.278000 

11.577000 

49.278000 

2 

-2.964000 

-24.165001 

-2.964000 

-24.165001 

-2.964000 

-24.165001 

1 

0.401300 

4.650300 

0.401300 

4.650100 

0.401100 

4.650  100 

4 

-0.028200 

-0.3141S0 

-0.023200 

-0.314180 

-0.028200 

-0.  114180 

'l 

0.524500 

0.524500 

0.524500 

0.524500 

0.524500 

0.524500 

6 

0. 640000 

0.640000 

0.640000 

0.640000 

0.640000 

0.640000 

7 

1 . 300000 

1.500000 

1 . 500000 

1 . 500000 

1 . 500000 

1 . 500000 

H 

225.899002 

225.899002 

225.899002 

225.899002 

225.899002 

225.899002 

9 

79.830002 

9.150000 

79.830002 

9.150000 

79.830002 

9.150000 

10 

1.570000 

1.570000 

1.570000 

1.570000 

1.570000 

1.570000 

11* 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

12 

12.000000 

12.000000 

12.000000 

12.000000 

12.000000 

12.000000 

n 

0.060000 

0.115000  ^ 

0.060000 

0.115000 

0.060000 

0.115000 

14 

6.000000 

6.000000 

6.000000 

6.000000 

6.000000 

6.000000 

15* 

0.000000 

0.000000 

0 . 000000 

0.000000 

0.000000 

0.000000 

15 

2.750000 

2.750000 

2.750000 

2.750000 

2.750000 

2.750000 

IT 

0.008789 

0.008789 

0.008789 

0.008789 

0.008789 

0,008789 

15 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

19 

0.004883 

0.004883 

0.004883 

0.004883 

0.004883 

0,004883 

20 

1 

1 

1 

> 

1 

1 

1 

21 

1 

1 

1 

1 

1 

1 

Parameter 

Tube  No. 

No. 

7 

4 

5 

10 

1  1 

12 

i 

H.  377000 

49.278000 

13.577000 

49.278000 

1 1.577000 

49.27800 

2 

-2.964000 

-24.165001 

-2.964000 

-24.165001 

-2.964000 

-24.165001 

i 

0,401300 

4.650300 

0.401300 

4.650300 

0.401300 

4.650300 

4 

-0.028200 

-0.314180 

-0.028200 

-0.314180 

-0.028200 

-0.314180 

5 

0.524500 

0.524500 

0.524500 

0.524500 

0.524500 

0.524500 

6 

0.640000 

0.640000 

0.640000 

0.640000 

0.640000 

0.640000 

7 

1.500000 

1.500000 

1.500000 

1.500000 

1 . 500000 

1.500000 

8 

225.999002 

225.899002 

225.899002 

225.899002 

225.899002 

225.899002 

9 

79.830002 

9.150000 

79.830002- 

9.150000 

79.830002 

9.150000 

1°  . 

1.570000 

1.570000 

1.570000 

1,570000 

1.570000 

1.570000 

11* 

0.000000 

0.000000 

0.000000 

0,000000 

0.000000 

0.000000 

1-' 

12.000000 

12.000000 

12.000000 

12.000000 

12.000000 

12.000000 

13 

0.060000 

0.060000 

0.060000 

0.060000 

0.060000 

0.060000 

14 

r.  000000 

6.115000 

6.000000 

6.115000 

6.000000 

6.115000 

15* 

: . 200000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

1  r» 

*52000 

2.750000 

2.750000 

2 . 750000 

2.750000 

2 . 750000 

Z .  IO?"  59 

3.09S789 

0.008789 

0.008789 

0.008789 

0.008789 

N 

3.  300000 

0.000000 

0.000000 

0.000000 

0.000001) 

0.000000 

1? 

j. 004833 

0.004883 

0.004883 

0.004883 

0.004881 

0,00488  1 

M 

1 

1 

l 

1 

1 

1 

1 

1 

l 

1 

1 

1 

*<*il#or.  Plot  values  not  entered 
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uncorrected 

1.035 


An  adjustment  for  axial  heat  loss  yields: 


h" 


h’ 


1  + 


2  I C  *  T 

lW  h' 


(5) 


(6) 


where  L  is  the  length  of  the  copper  block  in  inches,  C  is  thermal  con¬ 
ductivity  of  the  tube,  and  T  is  the  wall  thickness  of  the  tube. 

To  allow  for  comparison  between  values  of  h  calculated  at  different 
times,  h  was  referenced  tS^a  nominal  water  temperature  and  flow  velocity. 
The  nominal  water  temperature  and  flow  velocity  were  chosen  to  be  70°  F 
and  6  feet  per  second.  The  heat  transfer  coefficient  was  adjusted  with 
the  equations: 


...  *  1  +  0.0105*70 
"  1  +  Q.0105*T 


(7) 


where  T  was  average  water  temperature  in  degrees  Farenheit  measured 
during  the  current  test,  and 


£  -  pr  -  SLOPE  *  (F'°’8-6  °’8)  (S) 

where  SLOPE  was  the  slope  of  the  Wilson  Plot  line  and  F  was  the  average 
water  flow  velocity  measured  during  the  test. 

The  fouling  resistance  was  calculated  by  observing  the  difference  in 
the  rate  of  heat  transfer  from  that  of  a  clean  tube 


R  „  1 _ 1 

Rf  h  h 


initial 


where  was  t^e  initial  heat  transfer  coefficient  as  determined 

from  the  Wilson  Plot. 


Tube  data  analyses  were  performed  by  the  software  subroutines  HTCOEF, 
CCN,  and  LINFIT;  listings  of  these  programs  are  in  Appendix  A.  Figure  4 
shows  a  typical  analysis  printout;  Table  3  contains  an  explanation  of  the 
analysis  parameters.  These  parameters  are  also  stored  on  a  DEC  RKOl 
floppy  disk  for  future  retrieval.  Software  was  developed  to  print  out 
tube  data  in  a  more  easily  read  format,  and  to  compute  and  plot  the  daily 
average  and  standard  deviation  of  each  parameter.  An  example  printout  is 
shown  in  Table  4;  a  plot  of  mean  daily  fouling  resistance  values  is  shown 
in  Figure  5. 

(Text  Continued  on  Page  18) 
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TABLE  3 

ANALYSIS  PARAMETERS 
(Sheet  1  of  2) 


Parameter 

Parameter 

n 

Name 

Explanation 

i 

ITUBE 

tube  number 

2 

IDATE 

date 

3 

ITIME 

time 

4 

VZERO* 

minimum  voltage  to  which  thermopile  dropped  at 
the  end  of  the  cooling  curve 

5 

DEV 

standard  deviation  of  VZERO 

6 

IFIT, 11,12 

fit  number,  beginning  and  end  point  number  in 
cooling  curve  used  for  regression 

7 

VS TART* 

maximum  thermopile  voltage  at  beginning  of  cool¬ 
ing  curve 

8 

DEV 

standard  deviation  of  VSTART 

9 

TAU 

time  constant  (seconds)  of  cooling  curve 

10 

DEV 

standard  deviation  of  TAU 

11 

A 

intercept  of  line  in  logarithmic  cooling  curve 
plot 

12 

DEV 

standard  deviation  of  A 

13 

B 

slope  of  line  in  logarithmic  cooling  curve  plot 

14 

DEV 

standard  deviation  of  B 

13 

R 

correlation  coefficient  of  logarithmic  cooling 
curve  plot 

16 

RMS 

root  mean  square  error  of  cooling  curve  plot 

17-27 

same  as  6-16  except  for  second  fit  to  cooling 
curve  data 

2S 

T 

mean  water  temperature  (degrees  Celsius)  during 
cooling  cycle 

*Values  listed  in  A  to  D  converter  units  can  be  converted  to  voltages 
by  the  relation: 


NCSC  TM-271-79  ! 

Best  Available  Copy  j 

TABLE  3  ! 

(Sheet  2  of  2)  j 


Parameter 

Name 

29 

DEV 

standard  deviation  of  T 

30 

F 

mean  flow  velocity  in  feet  per  second  during 
cooling  cycle 

31 

DEV 

standard  deviation  of  F 

32 

HUNCOR 

uncorrected  heat  transfer  coefficient  (h) 

33 

HCOREM 

h  corrected  for  air  and  axial  heat  losses 

34 

TCOR 

temperature  correction  to  apply  to  h 

35 

HCREWT 

temperature  corrected  1/h 

36 

VCOR 

velocity  correction  to  apply  to  1/h 

37 

HNOM 

corrected,  final  h 

38 

FR 

fouling  resistance  (ft^-hr-°F  BTU**1) 

39 

IREC 

record  number  for  data  storage 

TUBE 

• 

DOTE 

601 

TIME 

i  i :  ,i  n 

V2£Ru 

283, 9 j 

DEV  2.53 

tit  i 

,'ROii  L'hTi 

F  G I  i  <  T  5 

VS TORT 

2237. 

IEv  2.8 

r.(U 

3o .  7  23 

SEv  0,107 

A 

7.7131 

DEV  0 . 000c. 

r* 

■■0,01497' 

DEV  0.00002 

R 

"0.97261 

RH5  17. 1 87 

r  it  2 

F'RGI’j  DrtT/j  FujlNT  6 

VSTAR*: 

2277. 

EEv  -j .  3 

Trt-J 

5  -1 .  6  -1  6 

DEV  0.223 

A 

7.7-10-1 

DEV  0.00*0 

B 

~0 , 01330 

DEV  0.0000-1 

R 

"1  .00001 

RMS  1.113: 

T 

28.271 

DEV  0.0 10 

r 

5.5-17 

DEV  0 . 086 

MUNCOR 

j  006 . 

HCOREU 

931. 

TCOR 

0 . 72.3 

HCKEwT 

835 . 

vCCR 

0,32327- 

0  1 

KMDN 

928. 

f'R 

0  .  j  E  0-1 

AR2H  WE 

r-  on  data 

OiSn  A£  ;.CvOR 

FIGURE  4.  TYPICAL  ANALYSIS  PARAMETERS  PRINTOUT 
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TABLE  4 

TUBE  DATA  LISTOUT 


CTRST  AVAILABLE  FREE  WCtK?  2056 
TUBE  NUMBER  2 


RE  El 

BATE 

TIM" 

CST 

y*KR0 
a.  r* 

TCMP 

C. 

EL  DU 

rps 

FL0WSP 

SEC 

TAU 

H.VnM 

R 

RMS 

RFOIIL 
E— > 

854 

622 

215. 

"'6  *?  * 

3.627 

0.031 

33 

.42 

707. 

- 1 . 00000 

. 7  6 

0 . 73 

958 

622 

30 

20". 

27.  13 

5.615 

0.023 

55 

•  76 

go  .*  ¥ 

-0.99994 

■1,63 

0.38 

962 

622 

i 

o 

203. 

24 . 98 

5 . 607 

0.003 

35 

.82 

39  4  . 

-0.99994 

4.9  4 

0.89 

866 

622 

i 

30 

206. 

*  21 

5.576 

0.023 

55 

.65 

901  . 

•0.99903 

4  . 59 

0.30 

870 

622 

*> 

0 

211  . 

27.12 

5.392 

0.003 

33 

.78 

896. 

-0*9999.9 

l  *  2  0 

0 . 86 

874 

622 

2 

30 

218. 

27  t  22 

3.678 

0,030 

55 

,48 

896 . 

-0.99995 

3.79 

0 .  £.3 

878 

622 

3 

0 

226. 

27« 

5.661 

0.000 

.36 

900. 

-i .00000 

7  *61 

0.81 

893 

622 

3 

30 

0*?A  # 

27.06 

5^648 

0.004 

53 

.30 

907. 

-0 . 99993 

3.18 

0.72 

SS9 

622 

4 

0 

20 7 . 

27.08 

5.603 

0.023 

35 

.46 

906. 

-0.99987 

3.51 

0.73 

e?: 

62" 

4 

30 

201  . 

27 . 02 

Zi  •  603 

0.012 

55 

.92 

891. 

-i ,00001 

3.86 

0.97 

8’ 6 

62? 

3 

0 

•~l  *>”7 

77.03 

3 » 5?5 

0.009 

55 

♦  j  7 

918. 

-0,99993 

2,62 

0,59 

901 

622 

5 

30 

232. 

27.05 

3.627 

0.009 

55 

.96 

886 . 

-0 ♦ 99990 

6  •  10 

0.98 

906 

622 

A 

0 

210. 

26.90 

3 . 636 

0.033 

rr 

.  19 

9 1 2 . 

**0 .  '’9997 

4.20 

0 . 67 

911 

622 

7 

10 

120. 

27.04 

3  *  77 A 

0.018 

54 

.96 

906. 

-0.99983 

6 . 96 

0 . 73 

9  i  o 

622 

7 

4ft 

95. 

27.  11 

5.790 

0.015 

54 

•  33 

918. 

-0,99982 

7.01 

0.59 

?27 

622 

8 

10 

90 . 

27 . 1 7 

r. .  709 

0 . 008 

54 

. 2.9 

927. 

-0 . 99986 

6 » 47 

0.48 

915 

62? 

8 

40 

139. 

27.16 

5.037 

0 . 026 

53 

943. 

-0.99994 

5 . 49 

0 . 30 

9 -’3 

622 

9 

10 

173. 

27.23 

5 . 807 

0.017 

54 

.26 

926. 

-0.99978 

6.60 

0 . 50 

9  51 

622 

9 

40 

182. 

27.22 

5.762 

0.003 

54 

.54 

92  i  . 

-0.99988 

6.34 

0 . 56 

959 

622 

10 

10 

184, 

26.95 

5.747 

0.004 

54 

.94 

911  , 

-0.99987 

8.20 

0.68 

967 

622 

10 

40 

184. 

26,94 

5.796 

0.016 

54 

.73 

913. 

-0.99984 

7.07 

0 . 65 

?75 

627 

u 

10 

187. 

26.95 

5.750 

0.004 

54 

•  67 

920. 

-0.99983 

7 .  1  7 

0 . 57 

9G3 

622 

11 

40 

189. 

27.15 

5.767 

0.003 

54 

.74 

9  14. 

-0 . 99984 

6.94 

0.64 

991 

622 

12 

10 

187. 

27.24 

5.796 

0.022 

54 

.52 

9.10. 

-0.? 9988 

7,5* 

0 . 59 

??? 

622 

1  2 

40 

185. 

27.45 

5 . 300 

0.019 

54 

.65 

910. 

-0.99986 

6.71 

0.6? 

.1  007 

62.; 

13 

10 

186. 

27.45 

3.743 

0.007 

54 

.48 

<?r%2  • 

-0.99992 

7.32 

0.54 

1015 

622 

13 

40 

192. 

27.43 

5.S09 

0.021 

54 

.53 

914. 

-0.99989 

6.97 

0.64 

ICO  3 

627 

i  4 

10 

192. 

27.58 

5.751 

0 . 003 

54 

.56 

919. 

-0.99981 

7.28 

0,59 

1031 

622 

14 

40 

193. 

27.74 

5.800 

0.005 

54 

.69 

906. 

-0.99982 

7.22 

0.73 

1039 

622 

15 

10 

194. 

27.97 

5.768 

0.025 

54 

.63 

9.0, 

•0.99983 

6.85 

0 . 69 

104  7 

622 

15 

10 

197. 

28.07 

5.752 

0.010 

54 

.44 

917. 

-0.99904 

6,94 

0 . 60 

1055 

622 

16 

10 

200. 

28.20 

3 . 737 

0.021 

54 

.43 

916. 

-0.99980 

6  »  63 

0.62 

1063 

622 

.16 

40 

203. 

28.39 

5 . 75? 

0.006 

54 

.  42 

914. 

-0.99983 

6.16 

0,64 

1071 

622 

17 

1ft 

210. 

28.65 

5.729 

0.037 

54 

.48 

912 . 

-0.99979 

5.96 

0.66 

107^ 

622 

17 

40 

218, 

29.03 

5.832 

0.018 

54 

.36 

903. 

-1.00003 

5.41 

0.77 

1087 

622 

18 

10 

224. 

29.25 

5.794 

0.021 

54 

.28 

90S. 

-0.99981 

6  ♦  06 

0.72 

1095 

622 

18 

40 

223 , 

29,30 

5.798 

0.010 

54 

.25 

908. 

-0.99992 

5.05 

0.72 

1103 

622 
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SYSTEM  RELIABILITY 


The  OTEC  computer  system  at  Pnaama  City  demonstrated  Rood  relia¬ 
bility  throughout  testing,  with  brief  and  infrequent  computer  "down'' 
times.  Computer  hardware  performed  well  despite  round-the-clock 
operation  and  an  operating  environment  (a  trailer  on  the  pier)  with 
poor  temperature  and  humidity  control .  Only  one  major  breakdown  was 
experienced  last  year  due  to  a  lightning  strike  blowing  out  one  of  the 
computer  power  supplies. 

Computer  software  was  continuously  modified  to  allow  for  system 
expansion  and  enhancement  of  control  techniques,  and  allowing  for  more 
accurate  analysis  of  data.  Additional  software  routines,  written  in 
BASIC  language  for  use  on  a  Hewlett-Packard  9830A  programmable  calcu¬ 
lator,  have  provided  for  cross-checking  of  results  as  well  as  allowing 
for  graphical  representation  of  data. 

Fouling  resistance  was  determined  for  each  tube  at  the  end  of  the 
cool  down  cycle;  a  maximum  of  48  analyses  were  done  for  each  tube  in  a 
24  hour  period.  The  main  loss  of  data  was  due  to  a  tube  being  dropped 
from  a  test  because  of  flow  or  heater  problems;  in  such  a  case,  no 
analyses  were  done  for  that  tube.  All  analysis  data  were  rechecked  and 
not  used  if  other  problems,  such  as  high  flow  deviation  or  high  RMS 
errors,  were  detected.  Errors  on  the  analysis  parameters  from  one  run 
are  listed  in  Table  5.  These  errors  varied  between  tubes  but  were  con¬ 
sistent  for  each  tube.  They  appeared  to  be  small  enough  to  allow  for 
adequate  day-to-day  comparison  of  RFOUL  values. 
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NCSC  TH- 271-79 
APPENDIX  A 

ANALYSIS  PROGRAM  LISTOUTS 


tt:«dxi:htcoef.for 

FILE  HTCOEF .FOR  3/22/79-2 
ASF  FOR  OTEC 


SUBROUTINE  HTCOEF< ITUBE > 

BYTE  ITUBE 

COMMON  /ADATA/I VTC ( 300 ) t I VTH (300 ) , IVFM ( 300 ) 
COMMON  /FLOPPY/INTSM)  »REALS(24> 

BIMENSION  Y(300) »UT (300) 


ITUBE’THE  TUBE  NUMBER 

COMMON  /PRE/I DTEST f IHR» IMINj NPRELM  > IEND t NZERO 
COMMON  /MISC/I5AMPI 
COMMON  /DATE/IDAY 

IDA Y» IHR » I MIN -DAY r  HOUR»  MINUTE  OF  START  OF  BATA 
N'LRELrf-tPOINTS  BEFORE  HEATER  OFF 
IEND- TOTAL  Nt  BATA  POINTS 
I SAMP-SAMPLING  INTERVAL  IN  SECONDS 

COMMON  /HPARAM/ACCN . BCCN , CCCN ? BCCN » RAB1 r RAB2 ,RAD3 r 

1  AKCYL . AKTUBE r HAIR i HINTR  f ALNGTH » TTHWLL  f  VNOH .  HRSLPE  r 

2  THINTf THSLFfFMINT.FMSl  F .IWN5T.IWN 


ACCi! t BCCN. CCCN t BCCN  CURVE  FIT  FOR  H=F (TAU) 

F.ntl  =  I/2  IB  OF  TUBE  (IN) 

RAB2J/2  OB  OF  TUBE  (IN) 

,c.i>I'3-i/2  OB  OF  TUBE  (IN) 

AK C >’ L-- THEE UAL  COMB.  OF  HEATER  CYLINDERS  ( BTU/HE  FT  F>. 
AKTUBE-- THERMAL  COND.  OF  TUBE  < BTU/HR  FT  F)  . 
HY»IR-C0H'7ECTIVE  COEFFICIENT  FOR  HEAT  LOSS  TO  AIR 
r-INTK- INTERC.  ON  ORDINATE  AXIS  IN’i/H  VR5.  i/V**.8  PLOT 
ALN3TH- LEN3TH  OF  HEATER  CYLINDER  SET  (IN) 

TTHL'LL-WALL  THICKNESS  OF  THIN-WALLED  SECTION  OF  TUBE  (IN) 
UNO, ‘•(--•NOMINAL  FLOW  VELOCITY  IN  TUBE  IN  FT/SEC 
H.RSLPE-' SLOPE  OF  l/H  TO  1/(V**.S>  CURVE  (WILSON  FLOT) 

THIHT -INTERCEPT  IN  CENT.  DEG.  OF  THE  THEF.MISTEF:  CURVE 
.  TWATER  =  TH:INT+THSLP*I  VTH ) 

.'H3LP-SL0FE  OF  THERMI3TER  CURVE 

FMINT •■INTERCEPT  IN  GAL/MIN  OF  THE  THERHI5TER  CURVE 
<  r LGW-FMINTf FMSLP*IVFM ) 

FMSLP--3L0PE  OF  FLOW  METER  CURVE 

IWN3T-SAMPLE  AT  WHICH  WINDOW  TO  START  AFTER  NPRELM 
TWN-NUMBEF:  OF  TIME  CONSTS.  FOR  FINAL  WINDOW 
RTHENT  FUDGE  FACTOR  TO  ALLOW  FOR  AIR  HEAT  LOSS 

7'UBEIO-ID  OF  HEAT  EXCHANGER  TUBE  (IN) 


COMMON  /ICOQL/IPRTC 


1 3  AMP  “  I  Si  iMP  I  /  50 

IF  TCDE  IS  TITANIUM  -  LOOK  AT  EVERY  OTHER  SAMPLE 
SO  ISAMP  SHOULD  BE  DOUBLED 

IF ( AKTUBE . LT . 50 . ) ISANP=ISAMP*2 
IMR~MC*D( IHR f  24  > 

INTS  (  i  )  --ITUBE 
IN~o(2)--IDAY 
INT3(3)=IHR 
INTS ( 4 ) -IMIN 

WRITE (7f 1300) ITUBEf IDAY» IHR fIMIN 
130:*  FORMAT ( /// '  TUBE  'fI2/'  DATE  'fI4/'  TIME  'fI2»':'»I2> 


A- 1 
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c 

C  GET  PARAMETERS  POR  THE  TUBE 

C 

TUBEID-RAD1*2. 

C 

C  GET  THERMOCOUPLE  ZERO  POINT  AND  STD. DEV. 

C 

VTCZRO-O. 

Il*IEND-NZERO+l 
DO  100  1*11 » IEND 
100  VTCZRO»VTCZRO+IVTC(I) 

VTCZRO-VTCZRO/NZERO 
STDZRO=0 . 

DO  110  1*11. IEND 

110  STDZRO*STDZRO+  < IVTC  < I ) -VTCZRO )  *#2 

STDZR0=SQRT(STDZR0/(NZER0-1 ) ) 

REALS< 1 )=VTCZRO 
REALS  <  2 )  “STDZRQ.^ 

WRITE  ( 7*  1000 ) VTCZRO  «  STDZRCt 
IF(STDZRQ.NE.O. ) GO  TO  606 
WRITE( 7.  6077) 

6077  FORMAT < // '  ANALYSIS  ERROR'//) 

RETURN 

606  CONTINUE 

1000  FORMAT (  '  VZERO  '.F6.2.'  DEO  '»F3.2> 

C 

C  SHIFT  THERMOCOUPLE  DATA  TO  THE  ZERO  POINT 

C 

NP:l»NPRELM+l 
DO  200  I*NF1 . IEND 

200  Y(I)*IVTC(I) -VTCZRO 
C 

C  PRINT  RAW  DATA  (THERMOPILE.  THERMISTOR.  FLOW)  IF  IPRTC=»1 
C  PRINTED  VALUES  ARE  IN  A/D  CONVERTER  UNITS 
C 

IF  (IPRTC.NE.l)  GO  TO  97B 
DO  201  1*1, IEND 

201  WRITE  (7.979)  IVTC ( I ) . IVTH( I > . IVFM ( I > 

979  FORMAT  (4X.3I12) 

978  CONTINUE 

C 
C 

C  SEl  UP  DATA  ARRAY  FOR  FIT 

C 

DO  300  I*NP1 . IEND 
ARG1*Y ( I ) 

ARG*ABS ( ARG1 ) 

IF  (ARG.LT, .001)ARG*.001 
YCI )*ALOG ( ARG ) 

IF  (ARG1  .LE.O.  )ARG1'*.001 
300  WT ( I >» ( ARG1/STDZR0  >  ##2 

C 

C  •  DO  FITS  OF  DATA 

C 

DO  999  IFIT*1 .2 
IF  ( IF IT. EG. 2) GOTO  350 
C 

C  FIRST  FIT  FROM  NPRELM  TO  IEND 

C 

I1*NP1 
I2*IEND 
GOTO  400 
C 

C  SECOND  FIT  FROM  IUNST  TO  IWN  TIME  CONSTANTS 

C 

330  I1»TWNST+NP1 
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12  ^11  + 1 WN*  TAIL  VSAMr 

c 

400  NPTS*I2-.Tt  H 

WRITE <7*273 l>IFir,Tl»l2 

2731  FORMAT <  T  IT  '  • 1 2 >  FROM  OAT A  POINT  '  *  TO  • 
l  '  TO  PAT A  POINT  ,13) 

CALL  LINFT T(Y< It ! ,WTtIl ) •  NPTS.A*SXO,'iAA*P* SIGMAS, 
1  P ,  RMS  « r.H  I  TOT ,  STD  ,  1  1 * .1 SANE,  NPKELM  r  ION  > 

C 

C  WE  HAP  VTC ( T  >" ( VO) ♦EXP(-T/TAU>  WE  HOPEP 

C  WE  HADE  V  <  r>»t.nO(VO'-T/1Ali 

C  OR  A  I  1HT 

r 

c 

c 

I'O-'-EXFirt) 

VOHtHH -E>TiAlStOHAA' 

volume  xr  t  a  -signaTP* 

TAIJ*-l  ,/P 

TAIjHI  -’-t  .  \  P  ISI'OMAP  • 

TAULOW-  -1./<P  OTOM AN ' 

INTAV  NPTl'.l  ISAHP 
PV-VOUTOK-VOi  OW 
D  T  AH  *  V  A  U  H  *  -  T  AH  l .  0  W 
REALS  !3 ''-•■VO 
REALS  v  4  >  *PV 

REALS  (3'  *TAU  >. 

REALS  <6>  -PI  All 
HEALS  <  7  )  r-A 
REALS  <  0  >  •■S.tdMAA 
REALS ( « ) - P 
REALS \  i0.'*8I0MAP 
REALS  <  I  l  >  --R 

WR I TE  <  7  *  .1 1 00 ' 00 .  PV ,  TAU ,  PTAU 
1100  rORMAT  V  •  VS  TART  •,F3.0*’  PEV  -,F3.i/ 

1  '  TAU  ••♦F7.3.  •  PEV  '»FA.3> 

WR  I  TE  <  7  •  1 1 30  >  A ,  3 IONA A .  P .  S 1  OMAR . K  r  RMS 
1130  FORMAT  (  '  A  '.F7.4*-  PEV  '  »r<S .  A.-1 

1  '  P  •»F0.S»*  PEV  '*77.5/ 

v-.  2  '  R  '  *  F8 » 3 ,  '  RMS  rFO.S> 

C 

C  OCT  WATER  TEMP  AND  FLOW  VCL ,  OVER  ANALYSIS 

c 

AVETW  ■  0  ,i' 

AVEFV  0, 

ALTNFM-.  40fl3/TUREIP**2 
PO  423  I«I1*I2 

AVETWn-AVETW+I  HINT LTHStFiT  V  ILK  I ) 
AVtFV-'AVEFV+ALINTM*<FMINT+FMSLF*IVFM<  t)  > 

•123  CONTINUE 

AVETW  =  AVETW  ,'K'PTS 
AVEFV*  AVEFV/Nf'TS 
8TPWT.-0. 
storv-o, 
no  450  i-:.t ,  is 

ST0WT'>8TP VT+  <  THTMT+THSLP*! VTK < I > -AVETW  >  * *? 

V  -ALIMFMt’IFMINT+FMSLPtTVFMI  I >  > 

-130  STDFV  ■  SI  PFV  I-  ( V-AVEFV  >  **2 

STPWT -'SCtST  <  STPKT /  (  NF  TS- 1  ) ) 

STPrV  SdRr(3TPFV/(NPTS-,l  J 

REALS! J 2)* AVETW 

REALS  < 13 > -STPWT 

POM  IM  I  4  ) ■"■AVEFV 

REALS  %  13) "  STOTV 

1 F  <  I  FIT.  FG .  :■  >  WR  tTE<7*  1 700 1  AVETW  *  STPWT  «  AVFFV .  STpr-'V 
1200  FORMAT <  ’  r  »,r7.3>  1  PEV  »F3,3-' 
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1  '  F  •  irA.i.  DEV 

c 

C  NOW  PET  H 

CALL  CCN<TAU*AVETW.AVEFV*Tr:T.RHS) 

IF  ( IFTT.  EQ.  2.'MRTTE(7*  1 300) 

1300  FORMAT  < /) 

C 

999  CONTI NilS 

C 

IFi .NOT. tOK) RETURN 
C 

C  WRITE  OUTPUT  TO  FLOPPY  DRIVE  IF  FOSS  IMS 

C 

COLL  ASSIGNOR*  ' OX  1 2 OU TP UT  .  DAT  '  ?  J  •  ‘ OLD ' 

DEFINE  FILE  P<2300*32*L‘r  IV) 

READ (6  I IIKECNO^* 
ircRrcMC’.or.osoo 
,i  ,or>;.  iKfrrw.LT.o) 

2  SO  TO  9099 

WRITE  <  8  '  IRECNO >  3.  NTS*  REALS 
WRITE (7*3077/ IRECNO 

3077  FORMAT (  '  ARCMIVFD  ON  DATA  DISK*  AS  RECORD  MiJMPER  '*!“*•'  /I 
IRECNO-'IRECNOH 
WRITE  ( 0 ' 1  >  IRECNO 

CALL  CLOSE  (13)  .  >• 

RETURN 

9999  WRITE!?* 2349) IRECNO 

23*19  FORMAT <  '  ILLEGAL.  RECORO  NUMBER  FOR  FLOPPY  OUTPUT  ’.nv 
X  .  XX*30<  XH*)// 

2  ‘  PREVIOUS  RECORD  NOT  ARCHIVED  ' 

3  /'  PLEASE  ARRANGE  FOR  NEW  DATA  DISK' 

<1  //IX  *30  (IK*)///) 

CALL  CLOSE! 8) 

RETURN 

RETURN 

END 

*  END 
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>r:  kyx.Pa! 

r  r l  ~  cc.'» .  f or  3/ r  i  ••  ?  •  -  z 
r.  asf  for'otd: 

C  MOtiFtCli  BY  RLS  ?  / 1 3  '  7ft 

w 

31)  V-  ROUT  INE  CCN  ( TAL1 » TWA  rs.~<  r  rL  WEL  .HIT.  r-:l  i S  ‘ 

C 

COMMON  /HPARAH/A  t  P » C  r  0  r  RADI  -  RAD2 » SAT-3  •  AKCYL  •  Ah TuBE  f 

1  HAIR  rFiXNTR  t  ALNGTIi  .  TTKWLL  .  VKOt!  . URS'.PH . 7K.TMT  .  TF51.P . 

2  Fill  FT  .  FKSL"  .•  IWK3T  .  I  UN 
COMMON  /FLOFPV/INTS:-t> » REALS i 24 > 

c  •* 

C  RAIR»0  DUE  TO  PREVIOUS  CORRECTION  IN  HTAUR 

C 

IF  I TAU.LT. ,0001) RETURN 
,  FFTA!j*ALOG<  TAU> 

~NH '  A4FN7AU*(P:hFNTAl'*(C4D*FNTAU>  > 

K'NCOR  -EXm-'FH ) 

C 

R"  HEh'T  --J  .  030 

HVETE-HUMCOR/RTHENT 

I*:.  .  /ALNC-TH*  J  2 .  )  ABORT ( AKTiiBElTTHWU  /'1 2 .  ) 

R WALLS  -WLLCON/SOR  T  i l-IVIjTC ) 

HCGRCW -HU3TI-:/  ( 1 . 0  f  RWALI.S  > 

AK'Dftlr  >,*1.4.0 10154  70 .  > /Cl.  I- ,  01  05*01. 8*  TWA  I  ER+32.>> 

HCREWT*- ANORN  f  HCPREW 
HCE'.vTF  *- 1  <  0/HCF.‘EWT 
C 

REALS  < 1 6  >  *KUNCOR 
REAL  S  <  j  7  >  -‘HCOREW 
REALSv IS) -ANORM 
REALS  a*7>-KCREWT 

r F  V I F IT .  EO  .  2 )  WR I  TE 0  7 . 100  > KUNCGR ,  MCOREto 
100  FORMAT <'  Hl'NCCJR  '.F3.0/'  KCOREK  SF5.0) 
r F  < IF I T . EO . 2  >  URI TE  0  7 , 200  > ANORM r  HCREWT 
200  FORMAT <'  TCOR  *»F6.3/'  HCREWT  '.F5.0> 

C 

V9R-FL0VEL** 0- . 8  > 

UN0N3R‘-VND,'i** < - . 0  > 

VELC0R-HR5LPE# <  V8R-VNuM8R ) 

KRNClM-HCEWTR-v’tLCOR 

HN0M«1./HRK0M 

RF0IJL*NRN0M-<HTNTR+HRSLFE*VNDK8R> 

REALS ( 20 J-VELCOR 
REALS*!  2.1  >~HNOM 
REAL SO 22 )~R FOUL 
REALS  C  23  >  •■'RMS 
REALSi  2-1  >*0.0 
r 

RF0ULP*RF0UL*i.E4 

IF(IFI7.E0.2)WRITE<7>  300  > VELCOR .HNOM  jRFUULP 
300  FORMAT*,-'  VCOR  '.Ell.  A/'  HNOM  ’.F5.0/'  FR  '  r  Fti ,  2.  '  E-04 '  >  . 

RETURN 

emb 
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TT:--»D<18LINrIT.=,a^ 

C  FILE-  LINFIT.FGR  6/15/7D-1 

C 

C  ASF  FOR  OTEC 

r 

SUBROUTINE  LINFIK VrWTfHPTS* Ar3IGMAArB» SIGMABrftr 
t  RM5rCMIT'JTr5TrirIJ  f  I  SAMP » NF RELH  t  IOK ) 

C 

DIMENSION  YU>»UT<1> 

C 

IOK- , TRUE . 

SUM  =0 . 

3UMX-0 . 

SUHY-O . 

SUMX2 -0 . 

SUHXY--0. 

SliMY2«0. 

I T I  ME  -  <  1 1  -2-NPRELM )  *  I  SAMP 

DO  50  I=i7NpTS 
X I  ■=  I  *  I S  AMP  t  I T I  ME 
YI--Y>I) 

WE  tGHT-WT  < I ? 

SUM "SUM+WcIGHT 
SUMX-SIJMX+ WEIGHT*  XI 
10000  FORMAT  < IX* 3E13.6) 

SUHY-"SUMY(  UEIGHT*YI 
SUMX2-SUMX2TWEIGHT*XI*XI  N 
SUMXY-SUMXY+WEIGHT*XI*YI 
SUHY2=SUMY2+UEIGHT*YI*YI 
30  CONTINUE 

c: 

E'El.TA-SUM*SUMX2-SUMX*SUMX 
A--  ( SUMX2*SUMY-SUMX*SUMXY ) /BELT A 
B~  ( SUMXY*SUM-SUMX*SUMY ) /DELTA 
C 

C  COMPUTE  RMS  ERROR  ON  COOLING  CURVE  FIT 
DELTAY--0 . 0 
DO  947  1=1 fNPTS 
XI”I*ISAMP+ITIME 
YI=EXP(B*XI+A> 

YJ=EXPCY<I) ) 

94 7  DELTAY---DELTAY+<YJ-YI>.**2 

RMS=5QRT (CELT AY/FLOAT <NPTS) ) 

r 

CHITOT =0 . 

DO  55  1=1 >  NPTS 

XXX= A +B*  Cl * I SAMP  + 1 T I ME ) 

CHITOT=CK I TOT+WT ( I ) * ( XXX-Y < I >  >  **2 
35  CONTINUE 

VAL 1  "SUMX2/I'EL7  A 
VAL2--SUM/DELTA 

999  IF(VAL1.GE.0.AND.VAL2.GE.0)G0  TO  606 

WRITE(7r6077> 

6077  FORMAT (//'ANALYSIS  ERROR'//) 

IOK= . FALSE • 

RETURN 

606  CONTINUE 

SIGMAA=SGRT ( VAL1 ) 

SIGMAB-SQRT (VAL2) 

VAL1 -DELTA* (SUM*SUMY2-SUMY*SUMY> 

IF(VALI .GT.O. )G0  TO  6061 
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WRITE(7r6077)  - 
I OK- .FALSE . 

RETURN 

6061  CONTINUE 

R= ( SUM*SUHXY-SUMX*SUMr  > /SORT ( WL i ) 

IF(SUM.EQ.O. )G0  TO  999 
IF  ( NPTS-2 . EQ . 0  >  GO  TO  999 
IF (CHI TOT . LT.O. . OR . SUM. LT . 0. 

1  .OR. FLOAT (NPTS) /(FLOAT ( NFTS )-2. ) . LT . 0 ) GO  TO  9990 

STD=SQRT ( (FLOAT( NF'TS )/ (FLOAT (NFTS) -2. ) )* ( CHITGT/SUM ) ) 

RETURN 

9995  WRITE ( 7 r  9996 )CHITOTr  SUM t NPTS 

9996  FORMAT < '  CHITOT  '.E14.6.'  SUM  '.E14.6.'  NFTS  '.13) 
IOK=. FALSE. 

RETURN 

END 

*  ^ 
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D 

2) 

3) 

4) 

5) 

6) 

7) 

8) 
9) 

10) 

ID 

12) 

13) 

14) 

15) 

16) 

17) 

18) 

19) 

20) 
21) 


VARIABLES  -  HTCOEF,  UNFIT,  CCN  Analysis  Programs  (others  defined  at 
beginning  of  HTCOEF  routine) 


VTCZRO  average  minimum  value  of  voltage  at  tail  of  cooling  curve 
(VZERO  on  analysis  print  out) 

STDZRO  standard  deviation  of  VTCZRO,  should  be  small 

NZEkO  #  points  used  for  VTCZRO  calculation 

NPRELM  #  points  at  beginning  of  cooling  curve  to  ignore 

NP1  first  good  cofiling  curve  point  (NPRELM  +1) 

IEND  last  cooling  curve  point 

11,12  first,  last  point  of  cooling  curve  to  put  through  linear  fit 

NPTS  if  points  put  into  linear  fit 

V0  beginning  voltage  (A/0  units)  for  linear  fit  (VSTART  on 

analysis  printout) 

DV  standard  deviation  of  V0 

TAU  time  constant  of  cooling  curve 

DTAU  standard  deviation  of  TAU 

A  intercept  of  logarithmic  cooling  curve  fit 

SIGMAA  standard  deviation  of  A 

B  slope  of  logarithmic  cooling  curve  fit 

SIGMAB  standard  deviation  of  B 

R  correlation  coefficient  of  linear  fit 

RMS  RMS  error  of  linear  fit 

AVETW,TWATER  average  water  temperature  (T  on  analysis  printout) 
AVEFV,FL0VEL  average  flow  velocity  (F  on  analysis  printout) 

I  FI  T  if  of  linear  fit  to  cooling  curve  (is  done  twice) 


22) 

ALINFM 

gal/min  to  ft/sec  conversion 

23) 

STDWT 

standard  deviation  of  water  temperature 

24) 

STDFV 

standard  deviation  fo  flow  velocity 

25) 

0UTPUT.DAT 

disk  storage  file  for  analysis  parameters 

26) 

IRECNO 

record  #  for  storage  on  disk 

27) 

SUM, SUMX, SUMY, SUMX2 ,$UMXY ,SUMY2  summation  indices 

28) 

ITIME 

beginning  sampling  time  minus  one  point 

29) 

XI 

time  for  sampT^  I 

30) 

YI 

point  #1  for  linear  fit  (normalized  logarithm  of  sample) 

31) 

WEIGHT 

weight  to  apply  to  YI 

32) 

FNTAU 

log  of  TAU  (base  e) 

33) 

FNH 

log  of  heat  transfer  coefficient,  L  (base  e) 

34) 

HUNCOR 

uncorrected  L 

35) 

RTHENT 

correction  to  L  for  air  heat  loss  inside  PVC  housing 
(horizontal  tubes) 

36) 

WLLCON,RWALLS  correction  to  L  for  axial  heat  loss  along  tube  wall 

37) 

HVDTE 

h,  corrected  for  heat  loss  to  air 

38) 

HCOREW 

h,  corrected  for  axial  heat  loss 

39) 

ANORM 

temperature  correction  for  h  (TCOR  on  analysis  printout) 

40) 

HCREWT 

h,  normalized  to  70%  temperature 

41) 

HCEWTR 

1/  HCREWT 

42) 

V8R 

average  flow  velocity  to  -.8  power 

43) 

VNOM 

nominal  flow  velocity  (6  fps) 

44) 

VN0M8R 

VNOM  to  -.8  power 

NCSC  TM- 27 1-79 


45)  VELCOR 

46)  HRNOM 

47)  HNOM 

48)  RFOUL 

49)  RFOULP 


velocity  correction  for  1/h  (VCOR  on  analysis  printout) 
1/h  (corrected) 
corrected  h 

fouling  coefficient  (FR  on  analysis  printout) 

RFOUL  *  1.0E4 
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